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BROWNIAN  APPROXIMATION  TO  FIRST  PASSAGE  PROBABILITIES 


I.  Introduction  and  Summary. 

Let  Xj,X2»..>  be  independent  and  identically  distributed 

with  mean  E(x.)  ■  y  .  Let  s  •  x,  ♦  ...  ♦  x  ,  and  for  a<0  <  b 
I  n  I  n  — 

define  the  stopping  times 

t  -  t(b)  -  inf(n:  s  >  b>  (t  -  t(0)) 
n  t 

and 

T  ■  T(a,b)  -  inf{n:  sn  $  [a,b]>  . 

The  probabilities 

(1)  P{T<m} 

and 

(2)  P(T  <  m,  sT  >  b} 

arise  in  a  variety  of  probability  models.  They  are  difficult  to  compute 
exactly,  but  under  certain  conditions  may  be  approximated  by  the 
corresponding  probabilities  for  the  Brownian  motion  process.  Siegmund 
(1979)  gave  a  heuristic  argument  based  on  Laplace  transforms  to  show 
that  this  approximation  can  be  improved  considerably  by  obtaining  what 
amounts  to  the  first  term  in  an  Edgeworth  type  expansion  of  (1)  and  (2). 
This  method  has  been  extended  by  Yuh  (1980)  for  studying  joint  prob¬ 
abilities  of  the  form 

(3)  P{t  <m,  sm  <  b-x} 

and 

(4)  P{T<m,  s_  >  b,  s  <  b-x) 

1  ■ 

as  well  as  related  conditional  probabilities,  e.g.  P(x < ■J*m* b-x)  , 


'  1 


which  arise  in  the  study  of  Kolmogorov-Smirnov  statistics. 

The  purpose  of  this  paper  is  to  give  a  direct  probabilistic 
calculation  of  a  one-term  Edgeworth  expansion  to  probabilities  like  (3). 
The  method  is  in  principle  applicable  to  (4)  although  the  computations 
are  much  more  involved,  and  no  details  are  given  in  this  case.  An  ex¬ 
ample  of  our  results  is  as  follows. 

2  3 

Theorem  1.  Suppose  y  -  0,  EXj  *  i,  and  y  *  Ex^  is  finite.  Let 
1  /2 

b  ■  Cm  .  If  the  distribution  of  xJ  is  strongly  non-lattice  in 

the  sense  that  lim  sup  |e  exp(itx  ) |  <  1  »  then  for  each  x  >  0  as 
|th- 

m 

p{t  <  m,  s  <  (C-x)m1^2}  *  I  -  ®(C+x) 

m 

(5) 

-  m”,/2<P(C-bO[28+  (y/6)(x2-?2-  1)]  +  o(nf,/2)  . 

Here  ft  -  E(s2  /2EsT  )  if  5  >  0  and  B  -  Es  if  C  *  0  ; 

<P  and  *  denote  the  standard  normal  density  and  distribution  functions. 

Remarks. 

(a)  Since  P(t  <  m)  -  P(s  >  b)  +  P(t  <  m,  s  <  b),  if 

—  m  m  — 

P{t  <  m,  sm  <  b-x)  were  known  exactly  for  all  x  >  0,  then  (at  least 
for  continuous  distributions)  one  would  obtain  P(x  ^  m)  by  letting 
x  0  .  Although  it  seems  plausible  that  (5)  should  hold  uniformly  in  x 
for  x  near  0  ,  and  in  fact  the  right  hand  side  of  (5)  with  x  *  0 
agrees  with  the  result  obtained  heuristically  by  Siegmund  (1979),  we 
have  been  unable  to  prove  this  uniformity. 

(b)  In  the  ease  c  *  0,  Theorem  1  is  equivalent  to  a  result  of 
Iglehart  (1974).  In  this  case  the  asymptotic  behavior  of  P{t+>m)  is 
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known  from  fluctuation  theory,  e.g.  Feller  (1966),  p.399,  which  shows 
that  (5)  is  true  with  x  ■  0,  although  a  completely  different  proof 
is  involved. 

(c)  Siegmund  (1979)  has  given  a  method  for  calculating  6  numerically. 

(d)  Under  the  stronger  assumptions  that  the  x's  have  a  finite  moment 
generating  function,  Borovkov  (1962)  gave  a  complete  asymptotic  expansion 
of  (3).  His  methods  use  complex  analysis,  and  the  results  are  not 
given  in  a  form  which  permits  simple  comparisons  with  (5).  Also 
Borovkov's  methods  appear  to  handle  the  case  x  ■  0  without  difficulty, 
although  they  do  not  seem  to  adapt  readily  to  two-sided  stopping  rules. 

(e)  Analogous  results  may  be  obtained  for  arithmetic  distributions. 

(f)  If  for  some  n  the  characteristic  function  of  s^  is  integrable,  _1_ 

one  can  obtain  a  similar  expansion  for  the  density  P{t<m,  sm  <  (£-x)m2}, 

which  can  be  formally  calculated  by  differentiating  (5).  As  an  applica¬ 
tion  one  can  improve  the  limiting  distribution  of  the  one  sample 
Kolmogorov-Smirnov  statistic  -  see  Yuh  (1980)  for  details. 

The  remainder  of  this  paper  is  arranged  as  follows.  Theorem  I  is 
proved  in  Section  2.  Section  3  discusses  the  case  EXj  *  0.  In 
Section  4,  in  the  much  simpler  context  of  Brownian  motion  we  describe 
an  alternative  approach  to  these  problems  and  use  it  to  rederive  the 
results  of  Anderson  (I960). 


2.  Proof  of  Theorem  1 

Let  F  denote  the  distribution  function  of  s  ,  n  *  0,1,...  . 
n  n 

It  is  easy  to  see  that 


J_  1 

P{t  <  m,  s  <  (c-x)m2}  -  /  P{s  <  (c-x)m^  |  t,  s  }  dP 

(6)  1  1  T 

-  /  F  ((c-x)m2-s)dP  f-  /  F  (-  x  m2  -  R  )  dP  , 

{T<m}  n  T  T  {i<m}  “  T  * 


3 


1  /2 

where  R  ■  a  -  Cm  .  Similarly 

IB  T 

l 

.2, 


P{s  >  (C+x)ni  }  *  P{t  <  m,  s  >  (c+x)in  } 

*  m  *  *  “  q  * 

(7)  I 

2 


/  Cl  -  Fm_  (xm*  -  R)  ]  dP  , 
{x<m}  m  T  ° 


80 


(8)  p{t  <  m,  am  <  (?-x)in  }  -  p{sm  >  (c+x)in  } 


X  J_ 

2  _  v  *  ,  2 


/  [  1  -  F  (xm  -  R)  -  F  (-  x  n»  -  R  )  ]  dP  . 

{rim} 


The  customary  Edgeworth  expansion  applies  to  the  first  term  on  the 
right  hand  side  of  (8) ;  hence  the  remainder  of  the  proof  is  a  detailed 
expansion  of  the  integrand  in  (8)  and  an  asymptotic  evaluation  of  the 
resulting  integral. 

To  carry  out  the  following  analysis  it  is  technically  useful  to 

modify  (8)  to  insure  that  in  the  integrand  m-T  is  not  too  small  and 

R  is  not  too  large.  Let  m,  ■  mfl  -  (log  m)  ).  A  consequence  of 
m  i 

Lemmas  2  and  3  below  is  that  for  some  e  ->-0 

m 


1 

2. 


P{x  <  m,  sffl  <  (c-x)n>  }  -  P{x  <  m} ,  Rffl  <  m  effl,  s^  <  (?-x) m  }  +  o(m  , 


and  hence  by  an  argument  similar  to  that  leading  to  (8) 


2 


(9) 


P{x  <  m,  sm  <  (c-x)m  }  -  P{sffl  >  (C+x)m  } 

x  x 

“  /  }  [1  -  F  (  xm2  -R)  -  F  (-  xm2  -  R)  ]  dP 

.  d  .  Z  x  ID  T  IB  BIT  m 

{x<m. ,R  <  m  c  } 

I  ID  ID 


According  to  Petrov  (1972)  VI. 3  Theorem  3 

I  I  -I  -I 

FR(xn2)  ■  #(x)  -  (y/6  n2  )  (x2-l)  <p(x)  +  (l  +  |xj3)  o(n  2  )  , 
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where  o ( •)  is  uniform  in  x  .  This  may  be  used  to  expand  the  integrand 

in  (9),  and  a  subsequent  expansion  of  the  normal  distribution  function 

1  /2 

♦  by  Taylor's  theorem  show  that  uniformly  on  {t  <  m. ,  R  <  m  e  } 

1  m  m 

the  integrand  in  (9)  equals 

-I  -I  I 

m2(l.-T/g>)  2  (p(x/(l  -  r/m)2)  {2Rm+  (y/3)  [x2/(l  -  r/m)  -  l]  } 

_I  -I 

ofm  2e  R  )  o(m  2)  . 

IB  IB 

Suppose  initially  that  ?> 0.  By  Doniker's  theorem  t/m  converges 

-3/2  1/2 

in  law  to  a  distribution  with  density  function  ct  <p(C/t  ) .  By 

1  /2 

the  renewal  theorem  R  ■  s  -  Cm  converges  in  law  to  a  distribu- 

IB  T 

tion  with  density  function  P{s  >  y}/E  s  ,  and  according  to 

Siegmund  (1975)  it  is  asymptotically  independent  of  r/m  .  Moreover, 

by  Lemma  1  below  {R  }  is  uniformly  integrable,  so  these  limits  may 

m 

be  taken  under  the  integral.  Hence  the  integral  in  (9)  is  asymptotically 

-I  I  -1  i  -2  I 

a2  /  C26  ♦  <Y/3){x2/(l-t)  -  2<p(x/(l-t>2)i;t  \{Ut2)  dt  , 

0 

which  after  some  calculus  yields  (5) . 

In  the  case  c  ■  0  the  argument  is  similar  but  much  simpler.  It 

is  obvious  that  t/m  ■*  0  in  probability  and  R  ■  s 

m  t+ 

2 

Lenna  1  ■  E  st  <  •  and  ^  >  uniformly  integrable. 

2 

Proof.  That  Es^  <  *  is  known  from  random  walk  theory,  e.g.  problem  6, 
p.232  of  Spitzer  tl976).  From  renewal  theory  it  is  known  that 

v  |  «y  2 

P{s  -  b  x)  -►  (Es  )  J0  P{s  >  y)dy  ;  and  since  Es  <  •  ,  renewal 
t  —  T+  t+  r+ 

2 

theory  also  yields  E(s  -b)  -►Es  /2  E  s  ,  which  proves  uniform 

T  T+  T+ 

integrability.  (The  uniform  integrability  may  alternatively  be  proved 
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from  first  principles  by  an  elaboration  of  the  indicated  idea  for 
.  2 

proving  Est  <  »  .  This  may  be  the  best  approach  to  use  in  Lemma  4 
below.)  + 

1  -i 

Lemma  2.  For  each  e>0,  PfR^  >  tin  }  ■  o(m  )  . 

Proof.  By  the  Markov  inequality 

i  -I 

Vp. 

“2 

which  is  o(m  )  by  Lemma  1. 

-2 

Lenina  3.  Let  nij  *  m(l  -  (log  m)  )  as  in  the  proof  of  Theorem  1. 

Then  as  m  • 

P{m,  <  r  <  m,  s  < 

1  1  m 

and 

P{m  <  t  <  m,  s  * 

J  —  IS 

Proof.  By  Lemma  2 

I  II 

P{m  <  t  <  m,  s  >  (c*x)m2}  *  P{m.  <  x  <  m,  R  <  ~xm2,  s  >  (e+x)m2  } 

»  m  1  m  z  m  * 

-i  1  _i 

+  o(m  2)  <  sup  P{s  >  -~xm2}  +  o(m  2) 

“m^nim,  m~n  2  } 

f  2  1  7 

<  sup  P{(m-n)  s  >  -r-xlog  m}  ♦  o(m  )  , 

“m^nlm  m"n  2  ^ 

*2 

which  is  easily  seen  to  be  o(m  )  by  Nagaev’s  (1965)  improvement  of 

the  Berry-Esseen  theorem  or  by  the  related  result  of  Petrov  quoted 

earlier.  A  similar  but  easier  argument  shows  that  the  first  probability 

•1/2 

in  Lemma  3  is  also  o(m  )  . 


(C-x)m2}  -  o(m  2) 

i  -I 

(C+x)m2}  •  o(m  2)  . 
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When  EXj  4  0  ,  results  analogous  to  Theorem  1  are  more  complicated 
technically.  Although  it  is  probably  possible  to  formulate  a  comprehens¬ 
ive  theorem,  it  would  be  extremely  cumbersome.  Therefore,  in  this  section 
we  briefly  discuss  two  important  special  cases:  (i)  when  the  distribution 
of  Xj  can  be  imbedded  in  an  exponential  family,  and  (ii)  when  EXj  is 
a  location  parameter.  Although  the  treatment  of  these  two  cases  is 
slightly  different,  it  should  be  apparent  that  modulo  certain  technical¬ 
ities  the  methods  can  be  applied  to  other  similar  problems. 

To  consider  briefly  the  simpler  case  of  an  exponential  family, 

suppose  that  the  distribution  of  x,  is  given  by  Ffl(dx)  * 

exp[0x  -  <>(8) ]FQ(dx) ,  where  Fq  is  a  strongly  non-lattice  distribution 

having  mean  0  and  variance  1  .  It  is  easily  seen  that  if/(0)  *  0, 

♦’(0)  ■  Enx,  ,  and  ^"(9)  *  var.(x,).  By  taking  F  to  have  mean  0  , 

0  1  0  1  o 

$  has  been  standardized  so  that  lii’(O)  »  0  ,  and  thus  sgn  Enx,  ■  sgn  0  . 

0  1 

For  0  <0  (0  >  0)  it  will  be  assumed  that  there  exists  a 

O  1 

0j  >  0  (0Q  <  0),  necessarily  unique  for  which  tK9Q)  *  4>(0|)  . 

The  basic  identity  (6)  remains  true  and  in  the  obvious  notation 
becomes 

(10)  pJt  <  m,  s  <  (c-x)rn2}  -  /  F  (-xm2  -  R  )  dP  . 

0  m  {T<m>  0»m  T  m  9 

However,  (7)  must  now  be  altered.  Assume  to  be  specific  that  00 

and  write  6j  for  9  .  It  is  easy  to  see  from  the  exponential  family 

structure  that  for  every  function  h  >  0  such  that  h  I,  ,  is 

—  lx*n} 

B(Xj,  ...»  xq)  measurable  for  all  n  , 

/  h  dP.  -  /  h  exp{-(0.-0  )s  }dPfl  , 

{t<«>  9o  {t<«}  1  °  T  91 


and  hence 


(0-0  )b  4-  (6-0  )b  i 

■  °  P  {s  >  (?+x)m2}  -  e  0  /  [l-F  (xm2-R  ) ]dPfl 

o  {T<m}  o*  o 


-  /  tl“F  (xm-R  )]exp[-(9  -6  )R]dP.  . 

{T<m}  e0*m  T  m  1  o  m  6j 

From  (10)  and  (II)  one  obtains  the  following  analogue  of  (8): 

j  (0,-0  )b  j 

Pft  {t  <  m,  <  (c-x)m  }  *  e  P6  ^Sm  >  ^c+x^m  1 


-(0,-0  )R 
/  {e  1  °  m 

{x<m} 


[1-F  (xm2-R  ) ]  -  F  (-xm2-R  ) }dP_ 

0  ,m-x  m  o,,m  x  m  a, 

o  l  i 


A  similar  identity  holds  for  P„  ft  <  m,  s  <  (£-x)m2  }. 

o  m 

o 

Modification  of  these  identities  and  subsequent  expansion  as  in 

-1/2 

the  proof  of  Theorem  1  shows  that  if  (0j“eo^  *  ,  then 


P0  {t  <  m,  sffl  <  (C-x)m  }  *  exp[2£U+Bm  2>  ]  [  1  -  4>(t+x+£)] 

1  -i  -1 

-  e2U  m  2  <p(C+x+£){2e  +  (y/6)  U(C-x)  +  x2  -  2  -  l]}  +  o(m  2) 


For  0q  the  result  is  formally  identical  provided  £  is  defined  as 
+  1/2 

-  m  (0,-0).  The  details  of  these  calculations  have  been  omitted. 
I  o 

(For  the  ideas  justifying  a  version  of  Lemma  1  in  this  context,  see 
Siegmund,  1979.) 

Since  the  right  hand  side  of  (12)  makes  sense  provided  only  that 
the  x's  have  a  finite  third  moment,  the  exponential  family  assumption 
appears  to  be  much  too  strong.  However,  the  likelihood  ratio  of  the 
exponential  family  plays  an  important  role  in  the  derivation  of  (12), 
and  avoiding  its  use  raises  some  additional  technical  problems. 
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To  minimize  the  number  of  unpleasant  technicalities  it  will  be 
assumed  that  9  is  a  location  parameter  (which  is  further  restricted 
below  to  be  non-negative).  Hence  let  Fq  denote  a  continuous  strongly* non¬ 
lattice  distribution  function  having  mean  0,  variance  I ,  and  finite 

third  moment  y  .  Let  F.(x)  =  F  (x-0)  ,  and  let  F.  be  the  n-fold 

o  o  8  ,n 

convolution  of  FQ  with  itself.  Let  P0  denote  the  probability 
measure  under  which  x]*x2***'  are  ^dependent  with  common  probability 
distribution  F^  .  Except  for  the  indicated  convergence  in  distribution, 
the  following  generalization  of  Lemma  1  may  be  proved  by  the  method 
suggested  in  the  parenthetical  remark  at  the  end  of  the  proof  of 
Lemma  I.  The  details  are  omitted. 


Lemma  4.  4s  8+0,  f  s_  dP„  -*•  E„s  for  all  0  <  \  <  2  . 

-  .'-.TOOT  — 

lT+«»}  +  + 

-1  /2 

Let  t,  >  0  and  8  =  £  m  for  some  fixed  £  _>  0  .  Then  the 

1  /2 

P0  distributions  of  (s^  -  ?m  )  converge  to  the  distribution  with 

density  function  (EQs  )  *  PQ{s  >  y)  and  have  uniformly  integrable 
first  moments. 

-i/2  .  .  . 

Assume  now  that  9  *  C  m  for  some  fixed  £  0  •  The  identity 

(10)  remains  true  in  the  present  context.  However,  instead  of  (11) 

consider 

\_  i 

(13)  /  exp(-2£y)P  {s  G  (5+dy)m2}  =  /  [  /  exp (-2£y)F  (m2dy-R  ) ]dP 

(x,°°)  {tSm}  (x , 00 ) 

which  in  the  exponential  family  model  is  actually  equivalent  to  (11). 

From  (10)  and  (13)  one  obtains  the  following  analogue  of  (8): 


Pfi{x  <  m,  sm  <  (c-x)ni  }  -  /  exp(-2£y)P0{sm  6  (£+dy)mZ  } 


(IA) 


<x,») 

/  [  /  exp(-2£y)F  (m^dy-R)-F  (-xm^-R)]dP  . 

{T«ra}  (x,«)  9’m  T  m  6’°  T  m  6 


9 


With  the  aid  of  Lemma  4,  this  identity  may  now  be  expa^-d  along 
established  lines.  (Expansion  of  the  inter”-^i  on  the  left  hand  side 
of  (13)  and  the  inner  integr-.j.  on  the  right  hand  side  is  facilitated 
by  integration  by  parts,  application  of  Petrov's  theorem,  and  integra¬ 
tion  back  by  parts,  which  has  the  formal  effect  of  expanding  F  as 

0  ,n 

if  it  had  a  density  which  had  the  appropriate  local  expansion.) 

In  this  case  the  resulting  asymptotic  expansion  is  as  m  -*■  °> 

I  -i 

ry  a 

P0{t  <  m,  <  (4-x)m  }  =  exp{2£[?+m  (B+2yC?/3) ] } [ I  -  #  (4+x+C)  ] 

-i  -I 

-  m  2  e2^  <p(c+x+O{20  +  (y/6)[(x+02-42-  1+4?2]}  +  o (m  2)  . 

4.  Brownian  motion 

Our  first  approach  to  studying  the  problems  of  this  paper  involved 
a  different  method,  which  unfortunately  seems  to  be  more  complicated 
and  to  require  stronger  assumptions  than  the  method  of  Sections  2  and  3. 
On  the  other  hand,  it  seems  conceptually  simpler  and  can  be  used 
formally  to  construct  the  answers  to  problems  where  none  is  known  in 
advance.  By  way  of  contrast,  one  must  know  what  the  dominant  term 
in  (5)  is  in  order  to  add  and  subtract  it  in  (6)  to  implement  the  proof 
of  Theorem  1 .  If  one  proceeds  formally  to  expand  (6) ,  he  gets  the  wrong 
answer. 

To  illustrate  our  original  method  in  a  very  simple  setting,  we 
apply  it  in  this  section  to  give  a  simple  derivation  of  the  principal 
probabilistic  result  of  Anderson  (1960).  The  argument  is  constructive 
in  the  sense  that  the  process  of  successive  reflection,  which  is 
usually  the  difficult  part  of  two-boundary  problems,  is  accomplished 
in  a  purely  mechanical  way  by  means  of  a  simple  recursion. 

Let  {x(t),  0  _<  t  <_  I }  denote  standard  Brownian  motion  and 
for  <  0  <  and  ^  +  ti  <_  g2  +  n2  def^ne  T  *  inf {t:  X(t)  = 

^i +  °iC  ^or  i  a  1  or  2  } .  We  shall  calculate 
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(15)  p{t  <  1,  X(T)  -  c2  +  n2T  I  X(l)  -  V } 

(except  for  the  case  p  «■  Sj  +  hj  ■  C2  +  ^2)  *  'rhere  *■*  no  ^osa  °£ 
generality  in  assuming  P  <  ?2  +  n2  ,  because  in  the  contrary  case 
one  can  use  the  argument  to  calculate  the  complementary  probability, 
namely 

P{T  <  1,  X(T)  -  ? j  +  n,T  |  X(l)  -  y}  . 

Let  P^  denote  the  conditional  distribution  of  X(t) ,  0  <  t  <  I 
given  that  X(l)  ■  p  .  Then  the  probabulity  in  (15)  equals 

Py{T  <  1,  X(T)  «  ?2  +  n2T}  . 

Let  F(t)  *  8(X(s),  s  t)  .  For  any  t  <  1  and  Pj  *  p  the  measures 

P  and  P  contracted  to  F(t)  are  mutually  absolutely  continuous 
y  yj 

with  an  easily  computed  likelihood  ratio: 
dP 

(16)  —  -  exp{(y-yj)[X(t)  -  jt  (y+Uj)]/(l“t)}  *  L(X(t),t;  p.Pj)  , 

U,,t 

say.  Hence  by  standard  likelihood  ratio  (or  martingale)  arguments 


Pjj{t  <  1,  X(T) 


C2*n2l} 


■  /  L(C-+n  T,T;p,p  )dP  . 

{T<l,X(T)-52+n2T}  *  yl 


From  (16)  one  easily  sees  that  the  choice  of  Pj  for  which 

i*2  +  n2T  ~  jT(y+Pj)  -  C2 ( 1  ~T)  ,  i.e.  Pj  -  2(c2  +  n2)  -  P  leads  to 


(17)  Pjj{t  <  1 »  X(T) 


?2+n2^"  exPf"  2  52(52+n2-p)]Pu  {t  <  I,  X(T) 


^2+n2 


T} 


Since  this  choice  of  p^  exceeds 
tion  p  <  ;2  +  n2)  ,  P^  {T  <  l}  - 


?2  +  h2  (by  virtue  of  the  assump- 
1  ;  and  hence  (17)  may  be  rewritten 
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(18)  Pu{T<l,  X(T)  =  C2+m2T}  -  exp[-2?2(52+n2-w)  ][I  -Pp  {t  <  1 ,  X(T)  -  c^T}!  . 

The  identity  (18)  may  now  be  used  recursively  to  calculate  (15),  since 
by  the  same  argument 


; 

I 

l 

\ 

> 

f 


(19)  P^  {t  <  1,  X(T)  .*>  C,+n,T}  >=  exp[-2c1(ci+n,+p-2(c2+n2))] 

[1  -  P^{T  <  1,  X(T)  =  C2  +  t,2T}]  , 

where  m2  =  2(ct  +  “  ?2  -  n2)  +  h  <  (Cj  +  hj)  .  etc. 

In  general  the  result  of  carrying  out  this  computation  is  an  infinite 
series  of  terms  which  alternate  in  sign.  For  the  very  special  case 
^  +  rij  *  C2  +  ^2  +  P  ,  it  turns  out  that  p2  ■  u  ;  and  it  is  only  necessary 
to  solve  the  two  equation  (18)  and  (19)  simultaneously  to  obtain 

exp{-2c  (C9+n5-p) }-! 

Pu(T  <  «T)  -  Vn2T)  -  ««pt2(t2-{,)(C2^2-U)')-l 

(cf.  equation  (4.24)  of  Anderson,  I960  for  this  result  as  well  as  the 
answer  in  the  general  case) . 
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